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Abstract 

Using an analytical variational approach we calculate the hole-hole contact interaction on the Neel 
background. Solution of the Bethe-Salpeter equation with this interaction gives bound states in d- and 
p- waves with binding energies close to those obtained by numerical methods. At t/J > 2 — 3 the bound 
state disappears. In conclusion we discuss the relation between short range and long range interactions and 
analogy with the problem of pion condensation in nuclear matter. 



1 Introduction 



Due to connection with high-Tj, superconductivity, the problem of mobile holes in the t-J model has received 
considerable development during recent years. This model with less than half-filling is defined by the Hamilto- 
nian 
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Where dj,^ is the creation operator of a hole with spin a (cr I) at site n of a two-dimensional square lattice. 
The operator dj^^ acts in the Hilbert space where there is no double electron occupancy. The spin operator is 
Sn = ^d\^(Tapdnj3. < nm > are the neighbor sites on the lattice. The Hamiltonian with a ~ 1 corresponds 
to the t-J model, and for a = it is the t-J^ model. Below we set J = 1. 

At half-filling (one hole per site) the t-J model is equivalent to the Heisenberg antiferromagnet model ||, ||] 
which has long range Neel order in the ground state ^ ^ . The problem is the behavior of a system under 
doping by additional holes. Strong interaction of mobile holes with localized spins makes this problem very 
complicated. 

Intense studies has lead to quite firm establishing of one particle properties in the t-J model. A single 
hole is an object with a complex structure due to virtual admixture of spin excitations. It has been shown 
both analytically and numerically [0-(l8) that one hole has a ground state with a momentum of 

k = (±7r/2, ±7r/2). The energy is almost degenerate along the line cosfc^: -I- cosky = which is the edge of the 
magnetic Brillouin zone. However, even the two hole problem remains controversial because there are many 
sophisticated effects in such systems. 

The hole can interact with spin waves. For the Neel background, the spin wave spectrum is gapless according 
to the Goldstone theorem. The effective coupling constant for the interaction of a hole with long wave length spin 
excitations for t/J<5 was calculated in our work Earlier it had been done in Ref.[^ in the perturbation 
theory Hmit {t/J <C 1) and for zS ^ 1 {z is the number of neighbor sites). 

If interaction of the Goldstone excitation with particles exceeds a certain critical value, the system becomes 
unstable. This was understood a long time ago for electron-phonon interaction (lattice instability) ||2l[. Just 
because of this fact and because the hole-magnon coupling constant is large enough, any small but finite doping 
destroys the long range antiferromagnetic order in the t-J model. It was shown for i <C J in Ref. For the 
Hubbard model, the instability has been proven in the Hartree-Fock approximation For the t-J model, the 
instability of long range antiferromagnetic order for t/J < 5 was demonstrated in our work |Q. In the same 
paper, we pointed out the direct analogy of this instability with pion condensation in nuclear matter (for review 
see Ref. p5|). 

In the present work we calculate contact interaction between two holes on the Neel background. "Contact" 
means that it is due to exchange by spin excitations with momentum q ^ n. There is no retardation in contact 
interaction. We consider this interaction for both parallel and opposite directions of the holes spins. We calculate 
the vertex function which provides the most general description of interaction. For the most interesting case 
of opposite spins, there is an attraction a,t t/ J < 2 — 3 which gives a very shallow bound state. Our result for 
binding energy agrees with that of Ref. [ p6| obtained by a numerical variational method in a restricted Hilbert 
space. It agrees also with recent result of Monte Carlo simulation However, there is a disagreement with 
results of exact diagonalizations on finite-size clusters ut 



|3|. Following Refs. ^, we suspect 
that the discrepancy is caused by finite-size effects in exact diagonalizations. Recent calculation |Q of the 
binding energy on larger clusters (up to 26 sites) and the scaling with system size indicated a critical value of 
J/t between 0.3 and 0.5, in agreement with our number. 

We base our study on results obtained in . The suggested trial wave function of a single hole was of the 
form -^kcr — ^kcr|0), where |0) is the background and /i]^^ is the creation operator of a dressed hole 
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Here A„ — +1 for the spin- up sublattice and An = — I for the spin-down sublattice, (5 is a unit vector corre- 
sponding to one step in the lattice, N is the total number of sites. The hole energy corresponding to ansatz (||) 
was found in Ref.jl^ by variational method 
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where ^ 

Jk = - {cos kx+ cos ky). (4) 



Expressions for the parameters eoj ^Oi 2;, and y in terms of ground state correlators are presented in Ref.|13 
With the normahzation of one hole per the lattice of N sites the coefficients in Eq. (0) are as follows 
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The parameters X, Y, u, and v are also expressed in terms of ground state correlators |13|] . In the present work 
we consider two backgrounds: 1) Ising background which is the ground state of the t-J^ Hamiltonian and has 
no quantum fluctuations {S^ = ±i). 2) Neel background which is the ground state of the Heisenberg model 
[a = 1 in the Hamiltonian (|l|)]. For these states the numerical values of parameters in Eqs. (^) and (H) arc 

I : eo = 1.0, Ao = 1.50, x = 0.00, y = 0.00, X = 1.0, Y = 1.00, u = 0.00, v = 0.00; 

N: eo = 1.2, Ao = 1.33, a; = 0.56, 2/ = 0.14, X = 0.8, r = 0.72, u = 0.42, w = 0.12. (6) 

Ansatz (||) includes states with a simple hole and states with one hopping to a nearest neighbor site {S^dl^_^g^ 

can be represented as —dl^_^_g^dnidl^i.). Using the string picture suggested in Ref.js^, one can say that ansatz 
includes only strings of the minimal length L = 1. It is quite natural for t ^ 1 , but it is a rather rough 
assumption for the Ising background at t > 1. However, for the Neel background this ansatz is justified by 
comparison with numerical simulations. It turns out that it provides good results up to t « 5 (see discussion 



in Refs. |g_3|, 18 ). It means that the hole motion on the Neel background is mostly determined by quantum 
fluctuations or, in the string language, that fluctuations cut strings. As far as we understand, a similar conclusion 
was made in Ref. [Q. Thus we believe that the wave function we use is valid for t < 5. 

The plan of the paper is as follows. In Sec. II we study the interaction of holes with parallel spins. In Sec. 
Ill we consider the case of the holes with opposite spins in the Ising limit and the bound state of two holes. An 
analysis of the Neel case is performed in Sec. IV. Sec. V presents our conclusions. 



2 Interaction of the holes with parallel spins 

We consider a state with two holes 

Using Eq. (|^), one can easily check that the effective creation operators h^^^^ obey the usual anticommutation 
relations: 

{4,., 4.4 = 0. (8) 

It means that in the present description the dressed holes are fermions and the two hole state (0) is well defined. 

In the naive perturbation theory, the vertex corresponding to scattering 1,2 3,4 is just the matrix 
element (3, 4|i/|l, 2). However in our case, the holes are composite objects and we must take into account that 
the overlapping (3, 4|1, 2) is not equal to zero. In this case the vertex is equal to 

r(3,4;l,2) = (3,4|i/|l,2)-(ei+e2)(3,4|l,2), (9) 

where e; is the single particle energy (^. Due to the energy conservation in scattering, ei + £2 — 63 + 64. This 
procedure with subtraction is well known in the many-body perturbation theory (see, e.g., Ref. ]3^). The 
physical meaning of Eq. (||) is especially evident at |3,4) — |1,2) when the interaction is just the difference 
between the total energy and the sum of single-particle energies. 

We consider now the Ising background (the t-J^ model). Within the accepted approximation, there is no 
dispersion for this background [see Eqs.(^) and (|6|)]. Therefore, we may introduce an effective creation operator 
/ijjg. of a hole localized at site n. Equations (|^|P'^reduce to 
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In the coordinate space the vertex has the form 

r||(TO', n'; TO, n) = (to', n'|i/|TO, n) — 2e{m' , n'\m, n), 



(11) 



where \m,n) = hl^^hl^^\0) (to ^ n and to and n are associated with the same sublattice). One can easily check 
that matrix elements (to', n'|iJ|TO, rt) and (TO',n'|TO, n) are not equal to zero only for to' = to and n' = n. 
(There is also a possibility of to' = n and n' ~ m which, due to anticommutation relation (^), gives the usual 
fermionic antisymmetrization.) Obviously, for distant sites to and n the second term in (11) exactly cancels out 
the first term. The interaction is not zero only for the closest possible to and n: (n, n + 2(5^; n, n + 25x) (Fig- 
la) and Fff (n, n + 5x + Sy\ n, n + Sx + Sy) (Fig. lb). The calculation is straightforward and leads to 



r|l(n,n + 2(5^;n,n + 24) = V|| = 4ti^^^(l - V^) - - 4// + 12// 



(12) 



Ti;]in,n + Sx + 6y;n,n + 6x+Sy) = = 16ij^/i^(l - 4/i^) - /i^ - 12/i* + 48/i^ 

The terms with t are due to the hopping part Ht of the Hamiltonian ([^) , and the other terms are due to the 
Hj part. In calculations we used that + Afj,'^ — 1 and e = v'^ + lO/i^ — Sti/fj,. 

The r||(TO, n; to, n) is in essence the potential energy in coordinate space. There is a very weak attraction 
at t < 1 which is due to configurations with closed holes minimizing (Hj) (see, e.g.. Fig. 2). In minima at 
t ^ 0.5, V|| w w —0.018. Practically within the accuracy of the calculation, it is equal to zero. At f > 1 the 
interaction (|l|) is repulsive. At t > 1, V[| w -^{t - ^) and Vj^ ^ j{t - |). 

One can easily convert the vertex into momentum representation [Eq. (|^) with |1,2) = 



r||(k3,k4 



ki,k2) = — 



V|| (COS^ Qx + COS^ Qy ) + V± COS qx COS Qy — {3 ^ 4) 
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Here q = ki — ka = k4 — k2; (^12.34 = 1 if ki + k2 = ka + k4 and (5i2,34 = otherwise. 

We expect that the vertex (13) calculated for the Ising background provides also a reasonable estimate for 
the case of the Neel background. Moreover, we think that for t > 1, due to the arguments presented after Eq. 
(^, it is actually better justified for the Neel background than for the Ising one. The calculation of the quantum 
fluctuation correction is rather cumbersome and we omit it for the case of parallel spins. We feel, having an 
experience with holes of opposite spins, that the correction is not large. 



3 Interaction of the holes with opposite spins on the Ising back- 
ground. Energy of the bound state 

Similarly to the previous case, it is convenient to use coordinate representation (|l0|). The contact vertex is 
defined by Eq. ( plj ) with the state \m,n) = hl^^hl^^\0) . For distant sites to and n, the second term in ( pi] ) 
exactly cancels out the first term and the interaction vanishes. As in the case of parallel spins, there is a 
potential-like term F||(n,n + 5x]n,n + 5x) (Fig. 3a). In addition, there appear effective hopping of clusters 
F||(n, n + n, + 5x) (Fig. 3b) and F|j(n, n — 5x]n,n ■\' 5x) (Fig. 3c) which turn out to be equal. Simple 
calculation gives 

Tu{n + 5x\ n,n + 6x) = U^ Atv^i{l - b^x^ + 4/*^) _ 1 _ 4^2 ^ _ ^2^6^ 

Ti^l(n,n - 5x;n,n + 5x) = F-fj^(n,n + Sy] n,n + Sx) = T, (14) 
T = 2ti^/i(l - V + 32/f4) _ ^2(^1 _ 4^2)(;l _ i2^x^). 

Again, the terms with t are due to the hopping part of the Hamiltonian ( [l|), and the other terms are due to the 
Hj part. To be accurate, we should say that within the wave function ( |10| ) in the t-Jz model there are more 
distant configurations which contribute to hole-hole interaction (Fig. 4). However, the corresponding values of 
interaction potential 

3 

r||(n, n + 2dx + Sy] n,n + 26x + 6y) = 3r-f^(n, n + 35x;n, n + '65x) = (1^) 



are very small and we neglect them. 

Now we can consider the bound state of two holes. It is quite obvious that in the present approximation the 
wave function of such a state with the total momentum p must be of the form 



jT^xi^n - r„,)hl^^hl^\0) exp[i-p ■ {vr, + r,n)]. (16) 



N 

The equation for bound state is given by variation of {ijjpKH — E)\ipp) with respect to the coefficients x(i)- If 
we separate the binding energy A {E — 2e + A), the equation becomes of the form 

5(^(^p|(i/-2e)|^p>-A(V'p|V^p)) =0. (17) 

The matrix elements (m', n'|(_ff — 2e)|m, n) are given by Eqs. (p^. In the term A(V'p|'0p) we set (m', n'\m, n) = 
^mm'^nn'- Actually, the correction to the last equality is small, and therefore it gives only a small correction 
to A. To avoid misunderstanding, we should note that in the effective interaction (12[l^) we calculated the 



overlapping {m' ,n'\m,n) exactly because 2e is much larger than the binding energy A. Further calculation is 
straightforward. Only the coefficients x(^) corresponding to the closest positions of clusters are not equal to 
zero. The obtained coefficients correspond to p- and d-wave states which turn out to be degenerate. The energy 
is minimal at p = 0, and the binding energy is equal to 

A = [/ - 2r = 16ti^fi^{l - 7fi^) - i - 2^2 _ ^ g^^e^ ^j^g^ 

Now we would like to derive Eq. ( p^ ) using the momentum representation. For the t-J^ model, it is just 
another mathematical way. Nevertheless, it is useful because for the t-J model, which we are mainly interested 
in, the coordinate representation is not convenient. Taking |1,2) = t^k-, i |Q) ^^'^ 1 3, 4) = ^t.'^\^t.4,i\^) 
Eq.(^ and using Eqs. (p^), one can easily transform the interaction ( p^ into momentum representation [cf. 
with Eq. (|13|)] 
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Here 7^ = 7ki. The simple structure of the vertex (19) reflects the simplicity of quasiparticle ansatz ([10[). The 



quasiparticles interact only when they are at the nearest neighbor sites which gives only 7q and 7^7^ terms. 
Let gk be the wave function of a pair with the total momentum p = 0: 



V'(r„,r„) = — iyj^gk/ik|/iLkilO>exp[ik- (r„ -r„)]. (20) 
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This function obeys the Bethe-Salpeter (BS) equation 

{E ~ 2ek)gk = rn(k, -k; k', -k')gk', (21) 

k' 

where the summation is carried out over the Brillouin zone (cos k'^ + cos k'y > 0) . For the Ising case — e 
[see Eq. (pO|)]. Using the vertex (19), one can easily find that there are two degenerate solutions of Eq. ( |2l| ) 
corresponding to two symmetries 

d-wave : gk = ■\/2(cos fc^, — cos /cj^), (22) 
p-wave : gk = 2 sin k^ ■ 

The B"fijj terms cancel in the process of integration. The value of binding energy is dertermined by the 7q 
part A — A = U — 2T which exactly agrees with Eq. (|l^) . 

Plot of the hole- hole binding energy A in the t-Jz model as a function of t is presented in Fig. 5. The 
result at i = (A = ~l/4) can be very easily understood. Two simple holes bound together have one less 
antiferromagnetic bond broken than when they are apart. As t increases, holes tend to part to gain delocalization 



energy. For t < 2 there is a reasonable agreement between our result and that of Ref . |2g] . However for larger 
t, our curve lies substantially lower than in Ref. . Keeping in mind that our wave function on the Ising 
background (|l^) is not very good for t > 1 [see discussion after Eq. (^], we would not like to insist on our 
result. Let us discuss the point in more detail. In our calculations for t — Jz model we make in essence an 
expansion in where z = 4 is the number of neighbor sites. Therefore, one can expect it to be valid up to 
t/ J ~ 4. Unfortunately at large t, there takes place a strong compensation in the hopping part of the binding 
energy ( p^ A — U — 2T ~ |(1 — j^. It means that this contribution is asymptotically 1/z. Corrections 

from including longer strings into the wave function ( |ic| ) are 1/z as well. Therefore they are very important for 
this value. The situation is different in the t-J model and now we come to this problem. 



Interaction of the holes with opposite spins on the Neel back- 
ground. Energy of the bound state 



Since we are finally interested in calculation of the binding energy, let us look at the BS equation (|21|). The Neel 
case differs from the Ising one in two ways: 1) There is a not-trivial dispersion [Eq.(|3|)]. 2) The vertex function 
r is renormalized by quantum fluctuations. Generally speaking, the calculation of fluctuation correction to T 
is straightforward: One should only substitute the wave functions (H) into Eq. (^. As we have pointed out in 
discussion after Eq. (|^), the wave function (^) on the Neel background is actually more reliable than on the 
Ising one. Unfortunately straightforward calculation of F is extremely cumbersome because there appear many 
complicated correlators. Thus we have to make some simplifications. 

1) For the Ising case the wave function of bound state (|2|) is concentrated near the edge of the Brillouin 
zone (7k — 0). We will prove that for the Neel background this concentration is even stronger. (The same 
conclusion can be made from the plots of wave functions presented in Ref. |^.) The reason of this effect is 
quite evident: The bottom of the single-hole band (||) is exactly at the edge of the Brillouin zone. Using this 
fact, we will consider the vertex F||(k3, k4; ki, k2) only for momenta lying at the edge of the Brillouin zone 
(7. = 0). _ ^ 

2) Another simplification is connected with calculation of background spin correlators 

{0\S^ S^S[ . . . |0). We make an expansion in the parameter a of the Hamiltonian (^ and restrict ourselves 
to the first order. Although the physical value is a — 1, this can be quite a reasonable way to estimate 
the influence of quantum fluctuations. For example, for the transverse nearest neighbor correlator this gives 
(0|S'+5^_^^|0) = 2qi = — ^ —0.167 which is in good agreement with exact numerical value. In the first order 
in a, it is the only correlator which differs from its Ising value. 

3) The last simplification is the following. We will calculate quantum fiuctuation correction only in the Ht 
part of the vertex. Actually, in the Ising limit there is no compensation in Hj contribution into the vertex ([l9| ) 
and into the binding energy (^8|). Therefore, one should to expect only small correction from fluctuations. It is 
not the case for the Ht contribution where happens a strong compensation: A = U — 2T ^ ^{1 — — > 
(see discussion in the end of the previous section) . 

Calculation with the given above simplifications is still cumbersome but straightforward. Using Eqs. (|^) 
and (H), we get the correction 



Aru(k3,k4;ki,k2) = -2git- 



lyfi•'il'iJc^ + 7ki+k3 + 7k2+k4) + (lei/"'// - 56j///)7, 



<5l2,34, (23) 



where 2qi = (0| 5+5,7+5-10) = — ^ is the nearest neighbor correlator. Clearly, the correction (g3|) is repulsive. 
The coefficients v and fj, in Eq. ( |23[) correspond to the Ising limit ( [lO| ) since the expansion in a is carried out 
explicitly. The i^fj!^ term in Eq. (^3[) comes from expansion of (3, 4|i?(|l, 2) in Eq. (^), the v^^^ term comes 
from overlapping (3,4|1,2), and the i/fi^ term comes from expansion of energy (ei +62). The (7ki+k3 +7k2+k4) 
part arises from interaction of the holes at next-nearest neighbor sites. Its relatively small value reflects the 
importance of quasiparticles interaction at the nearest neighbor sites. 

Adding the quantum fluctuation correction (^ to the Ising value of vertex (^, we get contact hole-hole 
vertex for the t-J model 



r'n(k3,k4;ki,k2) = — 
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A'^/d + S(7i73 + 7274) + — (7ki+k3 + 7k2+k4) 
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(24) 
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A' = A+ lOatiy^fi^, C = -atv^r-. 

3 



Expressions for A and B are presented in Eq. (|19|). In transformation from Eq. ( |23| ) to Eq. (|2J) we take into 
account the normalization condition i/^ -I- 4/i'^ = 1. 



Let us consider two hole bound state with zero momentum p = 0. If we expand the single hole dispersion 
near the band bottom 

1 2 At^lx + y) 

= + -A7k, ^ = [A2/4 + 4*^(1 + y)]i/2 ' ^^S) 

the BS equation can be rewritten in the form 

(A - A7^)5k = J2^'uiK -k; k', -k')gk', (26) 



where r'^ is given by Eq. (p4[). Let us first treat A-term perturbatively. At A = 0, in spite of appearance of 
(7ki+k3 + 7k2+k4) term in the vertex, the solutions (p2|) still satisfy Eq. (p6|) with binding energies 



A(,") ^ap = A'-C\ Af^ ^ad = A' + C. (27) 

The first corrections in perturbation theory are 

A^^^ = {9'P\>^llWp) - ^A, A« = {9f\\ll\9T) - ^A. (28) 

At A > 1 the corrections become too large and the perturbation theory is not valid. However, general solutions 
of Eq. (^6|) are obvious: 

sinfcj; 

p-wave: 5k oc (29) 

Ap - A7^ 

cos kx — cos ky 

d-wave : gk oc — — . 

Ad - A7^ 

Let us emphasize that the factor (A^ — A7^)^^ (s is the symmetry label, s ~ p,d) substantially enhances the 
wave function near the edge of the Brillouin zone (7k — 0). Substitution of the solution ( |2S| ) into the BS 



equation (26) gives the selfconsistency condition for A^. The binding energy vanishes at the point t = where 



Gs ( |27D vanishes. Near this point 



A,^-^^-{t-tsf (30) 



To be precise, we should note that for the d-wave there is an additional weak logarithmic (In dependence 
in the coefficient in this estimation. 

The numerical solutions for A^ as a function of t are presented in Fig. 6 (d-wave) and Fig. 7 (p-wave). 
Dashed lines correspond to a = 1 in vertex ( p4[ ) which is the physical value. There are higher orders in a which 
we do not take into account. To estimate their contribution, we present also results with a — 0.5 (solid lines). 
We believe that exact solution lies somewhere between a = 0.5 and a = 1. Of course the choice of a = 0.5 is 
rather arbitrary. However, we can recall the linear expansion in a for the parameters x and y which determine 
the single-hole dispersion (||): x — ~12qi = a, y ^ —Mi = "^/S To fit the physical values (^), we must 

take a « 0.5. Anyway, it is not very important since the interval between the curves with a = 0.5 and a = 1 in 
Figs. (6) and (7) is quite narrow. At most t, the d-state has lower energy. There is a good agreement between 
our results and results obtained by a numerical variational method |26) which are also presented in Figs. (6) 
and (7). Thus to our view, the bound state disappears at t w 2 — 3. This conclusion is in agreement with the 
results of recent works Q . 

Solving the BS equation for p 7^ 0, one can easily derive the dispersion of the bound state. For example for 
small p at t = J 

SEp{p) K0.l65p^, 6Ed{p)~0M8p^. (31) 

For p ^ the solution in general case has no definite symmetry (d or p) . The labels in above equation denote 
the symmetry of corresponding solution at p = 0. 



5 Conclusion 

We calculated the vertex function of contact interaction between two holes in the t-J model. Contact interaction 
is caused by exchange of spin excitations with momentum q ~ tt and has no retardation. We considered this 
interaction for both parallel and opposite directions of the holes spins. For the most interesting case of opposite 



spins, there was found an attraction at t/ J < 2 — 3 which gives a very shallow bound state (with p or d 
symmetry) . Our result for binding energy agrees with that of recent works p6| , E^, |3^ , but disagrees with 
that of exact diagonalizations on small-size clusters |jl^,|2^-|^. Following Refs. 0, ^6|, w e suspect that 
the discrepancy is caused by finite-size effects in exact diagonalizations. Recent calculation [ p4| of the binding 
energy for larger clusters (up to 26 sites) and the scaling with system size probably confirms this point of view. 

The real physical value oft/ J in high-T^ superconductors is J « 3 (see, e.g., ^ ^)- Therefore, the two 
hole short-range bound state obtained in the present work is probably irrelevant to high-T^ superconductivity. 
Moreover, we think that this bound state is by-product of the pure t-J model. Actually, the size of bound state 
is of the order of one lattice space a ^ las {as is the Bohr radius). Therefore, the Coulomb hole- hole repulsion 
on neighbor sites is ^ ~ e^/(ea) 4eV/e, where e is the dielectric constant. More realistic structure of hole 
wave function, which takes into account oxygen ions [^Sj, does not change this estimation. Even with the static 
value e ~ 50, the Coulomb repulsion F ^ 0.1 eV destroys the bound state ~ —0.3 J ~ —0.03 eV. Actually 
the Coulomb repulsion is even larger because in interaction we must use the dynamic value of the dielectric 
constant e(aj) with uo ^ Tc ^ 0.01 eV which is definitely smaller than the static value. 

Thus we have to introduce a short range hole-hole repulsion into the t-J model. One way to do it is to 
consider a t-J-V model |4^, However, we do not think that the results are sensitive to the specific way of 
introduction of short range repulsion. We believe that the main problem is long range dynamics at distances 
^ Pp^. This problem includes the long range antiferromagnetic order instability [ p4| . Nevertheless the short 
range repulsion is very important. 

The situation is quite similar to the pion condensation in nuclear matter p5[ , where both the interaction of a 
nucleon with the Goldstone excitation (pion) and short range nucleon-nucleon repulsion are important. At the 
moment we believe that suitable description of high-T^ superconductors is an effective long-range theory with 
spin-1/2 holes and Goldstone gapless spin-waves. (To avoid misunderstanding, we should note that long range 
instability destroys antiferromagnetic order and therefore Goldstone spin-waves definitely are not the physical 
excitations. They just provide a suitable basis set.) The single hole properties as well as hole - spin-wave 
vertex [|l^ are described by the t-J model. However besides that, we have to introduce "by hands" the contact 
hole-hole repulsion. 

It is interesting to notice that investigation of high-Tc superconductivity could help to better understanding 
of pion condensation. Actually the long-range instability (pion condensation) |25 and the BCS-type nucleon- 
nucleon pairing are usually considered as independent phenomena. However now, due to high-Tc supercon- 
ductivity, we understand that they are closely related. 
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FIGURE CAPTIONS 

FIG. 1. Contact interaction of the holes with parallel spins. Only the configurations with holes in the 
centers of square clusters are shown, but each hole hops within the corresponding cluster. 

FIG. 2. Configuration with closed holes which minimizes the Hj and gives a very small attraction 

between holes with parallel spins. 

FIG. 3. Matrix elements for the holes with opposite spins. Only the configurations with holes in the 
centers of square clusters are shown, but each hole hops within the corresponding cluster. Fig. a represents 
diagonal interaction, and Figs. b,c represent hopping of cluster. 

FIG. 4. Distant configurations for interaction of opposite spin holes in the t-J^ model. Only the 

configurations with holes in the centers of square clusters are shown, but each hole hops within the corresponding 
cluster. 

FIG. 5. Hole-hole binding energy in the t-Jz model as a function of t. The d- and p- waves are degenerate. 
Solid line is the result of the present work [Eq. (p^j. The circles are results from numerical variational method 
, and the triangles are results from exact diagonalizations [|T6[). 
FIG. 6. Hole-hole d-wave binding energy in the t-J model. Solid line is the result of the present work 
at a = 0.5. Dashed line is the same at a = 1. The circles represent the results of numerical variational method 




FIG. 7. Hole-hole p-wave binding energy in the t-J model. Solid line is the result of the present work 
at a = 0.5. Dashed line is the same at a = 1. The circles represent the results of numerical variational method 



